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Abstract: We explore the properties of non-local effective actions which include grav¬ 
itational couplings. Non-local functions originally defined in flat space can not be easily 
generalized to curved space. The problem is made worse by the calculational impossibility of 
providing closed form expressions in a general metric. The technique of covariant perturba¬ 
tion theory (CPT) has been pioneered by Vilkovisky, Barvinsky and collaborators whereby 
the effective action is displayed as an expansion in the generalized curvatures similar to 
the Schwinger-De Witt local expansion. We present an alternative procedure to construct 
the non-local action which we call non-linear completion. Our approach is in one-to-one 
correspondence with the more familiar diagrammatic expansion of the effective action. This 
technique moreover enables us to decide on the appropriate non-local action that generates 
the QED trace anomaly in 4 D. In particular we discuss carefully the curved space general¬ 
ization of lnD, and show that the anomaly requires both the anomalous logarithm as well 
as !/□ term where the latter is related to the Riegert anomaly action. 
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1 Introduction 

While the fundamental Lagrangians describing our known physical theories are all local, 
quantum loops of massless or nearly massless particles yield non-local effects. It is often use¬ 
ful to arrange those loop effects into a non-local effective action which enables a systematic 
investigation of the quantum effects on the classical background fields. For theories where 
the symmetries relate the couplings of different types of particles, such as chiral theories or 
general relativity, the evaluation of a single loop using the background field method allows 
the loop corrections to a large number of processes to be calculated at once. For example in 
chiral perturbation theory, the renormalized non-local effective action [1] is useful for many 
different reactions. 

In general relativity, Barvinsky, Vilkovisky and collaborators (hereafter referred to 
collectively as BY) have developed techniques for calculating and displaying the non-local 
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gravitational effective actions that arise due to graviton loops or those of other massless 
fields [2-11]. The results are presented using an expansion in the curvature. In effective 
field theory we are used to an expansion in the curvature for local Lagrangians. This 
corresponds to an energy or derivative expansion in which operators are suppressed by a 
mass scale which is typically the mass of the ’integrated out’ held. If the light fields present 
in the effective action are slowly varying, each term in the expansion is correspondingly 
smaller. Quantum mechanically, this corresponds to low energies. However, with non-local 
actions the curvature expansion has a different nature. Because non-local operators such 
as the inverse d’ Alembertian 1/V 2 can appear, higher powers of the curvature such as 
[(1/V 2 )i?] n are not automatically suppressed at low energy and the curvature expansion 
is not the same as the energy expansion. Instead, it is a way to describe the (calculable) 
infrared physics from quantum loops. The effects of these infrared non-local effects from 
loops are just starting to be explored [12-21]. 

In this paper we explore the non-local curvature expansion in a relatively simple setting 
- that of photons coupled to a massless charged scalar and to gravity. We also display the 
results relevant for massless fermions to highlight interesting features of the non-local action. 
Both the spacetime metric and the gauge field are treated as classical background fields. 
In a recent paper [22], we focused on obtaining the flat-space non-local effective action 
and the associated energy-momentum tensor that gives rise to the trace anomaly. Here 
we are concerned with generalizing the flat-space results to curved backgrounds. This is 
achieved via a technique that we refer to as the non-linear completion of the action where, 
similar to CPT, the action is displayed as an expansion in the curvatures. The non-local 
effective actions are a relatively unexplored topic and there remain interpretive issues that 
we explore here. 

Most notable is the issue of the covariant nature of the non-local form factors such as 
In V 2 . In particular, we pay special attention to the generalization of the flat d’ Alembertian 
to curved space which turns out to be a non-trivial aspect of the effective action. Moreover, 
direct use of the Feynman graph expansion of the effective action allows us to identify the 
terms which is related to the beta function of the theory and those which are not related 
to the latter. Our exploration leads to a better understanding of the non-local action that 
generates the QED trace (conformal) anomaly. To the best of our knowledge, this is an 
unsettled issue in the literature and the procedure of non-linear completion yields interesting 
insight into the correct form. 

The plan of the paper is the following. In section 2 we provide an overview of the 
main problem discussed in the paper and also present our results. In section 3 we discuss 
some of the methodological issues with this program, pointing out the main difficulties of 
constructing non-local actions in curved spaces and in section 4 we describe the non-linear 
completion matching technique. Section 5 is devoted for the non-linear completion of the 
quadratic action while the cubic action is displayed in section 6. We then move in section 
7 to show how the terms in the effective action generates the trace anomly. Finally, we 
summarize and conclude in section 8. 
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2 The problem of In □ 


In flat-space the one loop effective action for a photon, obtained by integrating out a 
massless charged scalar or fermion, has the form 


S 




bi In {U/p 2 ) 


F» v 


( 2 . 1 ) 


where bi is the leading coefficient of the beta function, b s = 1/(487t 2 ) for a charged scalar 
and bf = l/(127r 2 ) for a charged fermion, and □ = d 2 . Here the action is expressed in 
quasi-local form and the In a/p 2 operator is a shorthand for the fully non-local realization 


(x| In (H) |„) = *(.-„) = / In (=£) . (2.2, 


When one desires a formulation in curved spacetime, one requires that the logarithm 
generalizes to the covariant form, with tensor indices raised and lowered with the metric, 
and the □ operator also being covariant. We will reserve the notation □ for the flat-space 
d’Alembertian and use V 2 for the covariant version. That is, one requires 

j j d 4 x hr {a/p 2 ) F rx(i ^ j d A x^g g^g^F^ In {V 2 / p 2 ) F a p (2.3) 


This can be made more usable through the definition of the log as 


re 

In (V 2 //i 2 ) = — dm 2 

Jo 


V 2 + 


m * 


1 

p 2 + m 2 


(2.4) 


which then involves propagators that can be covariantly defined. Even here the result is 
not simple as the inverse operator is acting on the tensor indices of F a a and itself becomes 
a bitensor [23]. Later in the paper we expand the covariant form in eq. (2.3) to first order 
in the expansion = rj^ u + h^ Ul for photons satisfying p 2 = p' 2 , resulting in 


/ 


d^Xyfg F 01 ^ In 



F a /3 = 


d A x 


F a * In {a/p 2 ) F a0 + V (Of + Of) 


where 

Of = \vT v F^ In {U/p 2 )F af) - 2n log(D/ p 2 )F ua 

Of = dP&'F^F 0 * + d ll ff'F a3 ^F a s - rrd x F aS ^d x F a0 


(2.5) 


( 2 . 6 ) 


and indices are raised and lowered with the flat metric. We note that near the mass shell, 
p 2 = p’ 2 = A 2 ~ 0, the F{l/U)F terms are particularly dangerous as they involve the 
inverse photon “mass”. Notice also that the logarithms in eq. (2.6) are infrared singular. 

On the other hand, in our previous work [22], we have explicitly calculated the h fU/ 
corrections to the effective action for a conformally coupled scalar held and on-shell photons, 
and have extracted the fermionic analogy from the work of [24], Interestingly, none of the 
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above terms in eq. (2.5) are found in the result. Instead we get a relatively simple 
answer, in that the terms that are proportional to the beta function coefficient 1 are 


with 


and 


F a ? In ( n/fi 2 ) F a p + h 




2 In(D)r c 7 -~A IU ' 


= -F^FS + \ gia ,F af >F# 


A, 1U = d„F a pd„F a P + F a 0d lt d u F a P - ri^d x F a0 d x F a0 


(2.7) 


( 2 . 8 ) 


(2.9) 


This result is itself generally covariant to this order in h^, although different in structure 
from eq. (2.5). One can easily check that the full result is invarint under local coordinate 
transformations —> h^ + d(jAv) ■ I n contrast to eq. (2.5) we see that eq. (2.7) does not 

contain any of the dangerous F(1/\3)F terms - the inverse photon mass does not arise in 
perturbation theory. 

Both of the 0(h^ u ) terms in eq. (2.7) are required by trace anomaly considerations 
and hence must be proportional to the beta function coefficient The terms with the 
logarithm yield the correct trace anomaly for a pure scale transformation 

x' = Ax, A!^[x') = A _1 H M (x), h'^ y (x') = lnD' = lnD — InA 2 , (2.10) 

where the non-invariance of In □ leads to 

T£ = | {rr^F^p + 2 ( 2 . 11 ) 

which is the correct expansion of the covariant density ^fgF 2 . Under this rescaling the 
last term in eq. (2.7) is invariant. However, under a conformal transformation (g /iu —> 
exp (2a(x))g flv ) restricted to flat-space 

hfiu ^ hfiu + 2cr (2-12) 

the first two terms in eq. (2.7) are invariant while the last term is not. Using the on-shell 
condition \F\A^ = 0 we have that 2 

dxF a/3 d x F a P = ho (F ^), (2.13) 

and we see that last term yields the correct trace anomaly. The two related transformations, 
rescaling the coordinates and rescaling the metric, act differently in the effective action yet 
both yielding the same anomaly relation. We see that both types of non-locality, i.e. the 
logarithm and the massless pole in eq. (2.7) are required by direct calculation as well as by 
anomaly considerations. 

1 There is also a term independent of the beta function which we include below. 

2 The details of these steps are carefully presented in [22]. 
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We seek the covariant curvature expansion which reproduces the perturbative results. 
For nomenclature, the term of order F 2 is referred to as second order in the curvature, while 
that with an extra gravitational curvature, e.g. F 2 R , is called third order in the curvature. 
The details of the matching will be given in the body of the paper, while here we summarize 
the results. 

The mismatch of the two expressions eqs. (2.5) and (2.7) makes the expansion in the 
curvature relatively complicated. Because one is starting out with the Fln'V 2 F expression 
as the covariant form which is second order in the curvature, one needs to add and subtract 
correction terms in order to reproduce the actual calculated result. These counter-terms are 
third order in the curvarure as we show below. This does not modify the covariance of the 
result - both expressions are covariant. Nevertheless it does make the resulting expression 
at third order quite complicated. This matching procedue, which we refer to as non-linear 
completion , occupies most of the work described below. We find that the result to this 
order in the curvature is 3 


r log =i I d A x^\ F a p In (V 2 / fi 2 ) F a P - ^F afj F a ^R 


V 2 ' 


—I— 4 RL W — 

+ V 2 


log(V 2 ) t^-F^Ff + p^F a gF a A 


+ iVlog(V 2 )F/ - - A g^F afi \og(y 2 )F a P 


(2.14) 


where C a p is the Weyl tensor. Note that the logarithms within the square brackets [...] 
do not need a factor of p? as the log/x 2 would cancel between the two terms. In particular, 
these terms are scale invariant as we will discuss later on. 

We will show that eq. (2.14) has the correct anomaly properties. The way that this 
is accomplished is interesting. For a scale transformation as in eq. (2.10), it is the first 
term - the logarithm - which yields the anomaly. However for a local Weyl (conformal) 
transformation it is the second term - F 2 ^R - which is the active ingredient. This latter 
term appears as one of the portions of the Riegert anomaly action [25] when appropriately 
displayed in a curvature expansion. Finally for a global rescaling of the metric g^ u —> e 2rT g lliy . 
with a being a spacetime constant, there is a simpler path that again involves the logarithm. 
The latter is equivalent to a scale trasnformation as in eq. (2.10). We conclude that both 
the logarithm and the Riegert term (massless pole) are required by anomaly considerations. 
We comment on this dichotomy in regard to the geometric program to classify anomalies 
set forth by Deser and Schwimmer [26]. 

Finally in order to match the full result found in the direct one-loop calculation [22], 


li The placement of the differential operators appears somewhat different than the expressions in the 
body of the paper. This is allowed indeed under integration by parts as we are assuming asymptotically 
flat spacetimes. 
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one must add a nonanomalous term that has no relation to the beta function 


r weyi = nb I d 4 x^F^Fj>^C%„ 


(2.15) 


where 


n C = 


1 

967T 2 ’ 



1 

48-7T 2 


(2.16) 


This is different for fermions and scalars and is invariant under both scaling and conformal 
transformations. 

Our final result for the covariant one-loop effective action is 


where 


is the classical action. 


tot 


— SqI r log + T'weyl 



1 

Ae 2 (y) 


F^ v 


(2.17) 

(2.18) 


3 Covariant non-local actions: General remarks 

General relativistic actions are readily described when local. Using the metric, covariant 
derivatives and curvature tensors one can construct generally covariant local functions of 
the field variables. The ultraviolet divergences of quantum loops are therefore simple to 
treat because they are also local [27, 28]. However non-local objects are difficult to describe 
in a generally covariant form because they sample the metric at a continuum of points in 
spacetime. For a general metric, explicit expressions for such actions are not possible. 

For massless scalar QED and after integrating out the charged scalars at one loop the 
effective action must be gauge invariant and thus involves only the field strength tensor. 
Up to quadratic order in the gauge field and using dimensional regularization, a general 
form in curved spacetime is 

F[g,A] = \s E m+ [ d A x f d 4 y F^(x)M^(x, y, y)F a ^{y) (3.1) 

e o J J 

where Sem is the classical Maxwell action, eo is the bare electric charge and M^(x,y, y) 
is an antisymmetric second-rank bi-tensor density of unit weight which explicitly depends 
on the renormalization scale. As we show below, this bi-tensor samples the full space-time 
and not just the pair of points (x, y ) since it involves the effects of massless propagators. 
The practical question is what the form of this bi-tensor is and how we can best describe 
it. 

The divergence contained in M^(x,y) is local and calculable. It has the form [22] 

= 1927^6 ~ Vi v) ( 3 - 2 ) 
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where 



This divergence is absorbed into the renormalization of the electric charge. After removing 
this divergence, the residual bi-tensor L^i(x,y, //) is finite. 

One might expect that there are also local terms proportional to the geometric curva¬ 
tures, such as FFR which would correspond to M^ u ’ a d ~ <5 4 (x — y) R. Such terms are 

found when one integrates out a massive charged particle [29] . However, they are absent 
in our problem, that of integrating out a massless field, simply on dimensional grounds. 
The curvatures involve two derivatives of the metric, and hence the coefficient of any local 
term of the form FFR must have dimensions of 1/mass 2 . Because all fields are massless, 
there is no way to obtain such a coefficient. Any factors of the curvature in the action must 
be balanced by non-local factors such as 1/V 2 . This tells us that once we have dealt with 
charge renormalization, which is of course a local operator, the remainder of the effective 
action will be purely non-local. 

In flat space, the non-local function was obtained in [22] 

l( ° 17( x > V> M) = b ^~ I al3 L ( x ~ V’ d) (3-4) 

where b s = l/(487r 2 ) is the leading coefficient of the QED beta function for a charged scalar, 
e is the physical charge and L(x — y; //) is displayed in eq. (2.2). As a warm-up for later 
usage, let us pause at this stage to show how one can convert from a non-local form to a 
quasi-local one employing non-local form factors. The latter are the building blocks of the 
curvature expansion. Through the position-space representation 

(x| in (^) |y) = L(x - y) (3.5) 


one can re-write eq. (3.1) in quasi-local form as 


r ( 0 ) [A] = Sem + — 


d A x F^ u 





(3.6) 


To appreciate the subtleties in the construction of the bi-tensor, let us quote the effective 
action linear in metric perturbation around flat space g^ v = + h In non-local form, 

it reads [22] 


T«[A,/i] = - i I d*x [ d A yh^{x) 


b s L(x - y)Tji(y) - iy A F (x - y)T^(y) 


(3.7) 


where photons are taken to be on-shell, i.e. dropping factors of \FF, W . Here we have defined 


and 


= 2d ll F a pd u F°‘P - v^dxF^F^ 

T% = -F^FJ* + 


(3.8) 
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is the classical energy-momentum tensor. We also have the massless propagator 1 



(3.10) 


In this case, the result is local in the relative position of the gauge fields, but contains both 
logarithmic and massless-pole non-localities with respect to the gravitational field. Allowing 
the gauge fields to be off-shell would lead to a non-locality in all three held variables due to 
the appearence of the triangle diagram 5 . Let us arrange the bi-tensor density at this order 
in metric perturbation, it reads 



/ 


d 4 zh aX (z) L(z-xin)j£fa X + A F (z-x)H^ aX 5 {A \x-y) (3.11) 


where 



(3.12) 


One immediately notices that the bi-tensor density samples the gravitational held over the 


whole spacetime manifold. This is the main reason that the explicit construction of such 
non-local objects is not possible in arbitrary geometries. Instead, one can use the quasi-local 
form factors to express the loop correction as follows 



(3.13) 


where the position-space representation of the inverse d’ Alembertian is given above. In 


this paper, we seek a generally covariant non-linear completion of the above results that is 


accomplished by employing the non-local form factors. 

4 Non-linear completion: Expansion in the curvature 

The curvature expansion is a covariant method to display the effective action with arbitrary 
background helds. For local actions, the heat kernel expansion is the most elegant technique 
to resolve the functional determinant of any operator [5, 30-34]. Its usage encompases many 
applications in physics and mathematics, but unfortunately it becomes somewhat compli¬ 
cated when we deal with a massless operator. Moreover, the correspondance with the more 
familiar perturbative expansion of the effective action in terms of Feynman graphs is not 
very obvious [35, 36]. In this paper, we propose a new technique to obtain the effective 
action which we call non-linear completion. The logic is very similar to the matching pro¬ 
cedure well known in effective field theory (EFT). This procedure proceeds by perturbative 
matching of the full theory onto the effective theory. What makes the construction of the 

4 The boundary condition imposed on the propagator depends on the application one is considering. For 
instance, for time-dependent systems one should choose the retarted propagator. 

5 See the discussion in [22]. 
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EFT Lagrangian possible is the fact that it must inherent all the exact symmetries of the 
full theory. This is the pathway we are going to employ in our case as well. 

In our example, the symmetries of the full theory are diffeomorphsim and gauge invari¬ 
ances and hence the non-local action must be constructed from the generalized curvartures. 
As we have shown in the previous section, the form factors are an important tool as they 
enable the action to be written in quasi-local form where the action is manifestly covariant. 
One starts by listing the relevant curvature basis and organize it in terms of a power series. 
For the example at hand, we have 

1Z 2 : 

n :i : F^F^R, F m F v a W\ F^F^R^p, V^ v V a F v a . (4.1) 


The field strength is the curvature of the gauge-connection and thus counts as one power 
of the curvature. The effective action will be displayed as an expansion in these generalized 
curvatures. The last operator in eq. (4.1) does not contribute when the photons are on-shell 
and thus we are not going to discuss it further. Then one proposes all possible non-local 
functionals of the d’ Alembertian which could possibly act on the different terms in the 
curvarure basis 

/V 2 

* 1 ln (u 

p . J_ ] 

' v 2 ’ 

where the subscripts denote the curvature upon which the operator acts. As far as J -3 is 
concerned, one can arrange more operators such as 



ln(V 2 /V 2 ) 

/(V 2 ) 


(4.3) 


where /(V 2 ) is some function to be determined. However, we will see that no from factor 
of this kind arises in our example due to the on-shell condition. Although the above form 
factors look very complicated, these are all well defined via their Laplace transform 

/*00 

^(V 2 )= / dsF{s)e~ sV2 . (4.4) 

J 0 

The last step is perturbatively matching the full theory diagrams onto the non-local action. 
The ’Wilson’ coefficients in this case only depends on the coupling constants of the full 
theory and are to be adjusted via the matching procedure. Since a massless field is being 
integrated out, these coefficients can not depend on any mass or renormalization scale, i.e. 
the non-local action is completely insensitive to the UV. 


5 The 7Z 2 action: The elusive logarithm 

In this section, we discuss the non-linear completion of the flat-space action in eq. (3.6). It 
reads 

r [fh A ] = J J ^xy/gg^g^ F^ log (jj'j 
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(5.1) 




where V 2 = g^V^V v is the covariant d’ Alembertian. The matching onto eq. (3.6) is 
immediate. Now one must raise the question: what is the expansion of the above action 
around flat space? In partiuclar, the piece linear in the metric perturbation and its connec¬ 
tion to the perturbative computation. The answer to these questions is very important in 
understanding the covariant nature of the quasi-local expansion. In the remainder of this 
section, we show how to consistently expand the logarthim and prove that the 0(h) term 
in the action is entirely absent from the perturbative computation. We start by showing 
the steps for a scalar field as a toy example and then discuss the more interesting example 
of a 2-form. 

5.1 Toy example: A scalar field 

Let us consider the following action 

r b,0] = y d*Xy/g(j)\n 0 ■ (5-2) 

The goal is to expand the action around flat space to linear order in the metric perturbation 
+ h^ u . The most convenient way to accomplish this is to first vary the action 
with respect to the metric and then restrict the result to flat space. Using eq. (2.4), we 
find 

S g T[r],4>] = j d A x J dm 2 {</> (□ + m 2 ) 1 [5gS7 2 } g=ri (□ + m 2 ) 1 (j)} + ... . (5.3) 

where the ellipses denote terms resulting from the variation of yjg which do not matter to 
our discussion. To arrive at the above expression, we have used the formal variation of an 
inverse operator 

5g V 2 + m 2 = “ V 2 + m 2 V 2 + m 2 ' 

The variation of the d’ Alembertian depends on the tensor field in the action. For a scalar 
field, we have 

(^V 2 )* = (S^d^du - Sg^T^da - grST%,d a ) T (5.5) 

where 

= J g af} {d^5gp„ + d v 5gp^ - dgdg^) . (5.6) 

It is advisable at this stage to express eq. (5.3) in a non-local form which is accompliahed 
via the identity 

n j = / d 4 y A(x - y)^{y) (5.7) 

U + m z J 

where 

(□ + m 2 )A(x — y) = d^(x — y) . (5.8) 
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If we recall that Sg liu 


h around flat space, we find 


5 g T[rj,(j)] = J drxdrydrz J dm 2 4>(x) A(x — y) 

-h^dgdu - d ,l h, iv d 1 ' + ^ d a hd a ') A (y - z)(j){z) 


(5.9) 


where 


A(x - y) 


dH e —«■(*—») 

(27r) 4 l 2 — m 2 


(5.10) 


Although the above must be defined with some boundary condition, this is not going to 
affect our discussion. Notice that one could obtain the same result using the more explicit 
variation of the propagator 


6G(x,x') 

5g^ v {z) 



' £V 2 ' 
6g^ u (z) 


G{y, x') 


(5.11) 


To facilitate comparison with the perturbative calculation, we can Fourier transform the 
above expression and find 


f d A p d A p' 

J (2tt) 4 (2tt) 4 


<Kp)0(p / )^ /M '(<2 , )7V 


lnp /2 — lnp 2 

p2 _ p/2 ’ 


q = —p — p 


(5.12) 


where 


'PflV — 


1 , , 1 1 , 1,1 1 

2 Pflu + + 4 q i jPl ' + 1 q,/p i‘ + 1 q ^ Pv + 4 q,yP i J 


■ p'rinv 


1 

p ■ Whw 


(5.13) 


5.2 2-forms 


We now turn to the treatment of 2-forms which is our main interest. There are two distinct 
pieces that arise from the variation procedure. The first comes from varying the explicit 
factors of the metric tensor in eq. (5.1) while the second comes from varying the logarithm 
and the procedure is almost identical to the scalar example aside from some differences 
related to the tensor rank that we now discuss. First, we generalize eq. (5.5) to the 
variation of the d’ Alembertian when it acts on a 2-form 


(5 g V 2 A gu )\ g=ri = (~h ap d a dp - A gu - d p {ST^A au + 5T°p v A ag ) 

- 5Y^d p A au - 5T^d p A ga . (5.14) 

Second, we need to generalize eq. (5.7) 

□ + m 2 ^V = J A A g(z - y)Aap(y), - v) • (5.15) 

We recognize in eq. (5.14) a structure identical to the scalar held and the result is the same 
as before but with the difference that both transversality and on-shellness are taken into 
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account as described in eq. (A.7). We now show how to treat the new structures in eq. 
(5.14). In position-space, we have the following piece 


drxdrydrz / dm z A flu (x)A(x — y) 


~ 2 <^ r A A1 (y)(c )A A(y - x))A av {z) 


- (d x 6T^(y))A(y - z)A av {z) 


(5.16) 


where we used eq. (5.15). We now have all ingredients and after a laborious computation 
in momentum-space one finds 


(2) hi' d 4 u d 4 r/ 

r[<?,u] = r<»>M + jj (2a -)4 { ^ )i A °(p)A ll (-p')h , “'(-<i) -JW) + o(h 2 ) 

(5.17) 


where 

= \{q 2 r]^ ~ q^q„ - Q^Qv){p ■ PVa/3 ~ p' a Pp)—7T~ (5-18) 

z P P 

N^vaf} = M% a p log (5. H.9) 

and -M? a p is the tensor is the lowest-order matrix element describing the local coupling 
of photons to gravity. Explicitly, it reads 

= P'flPvVaP + PnP'vVa/3 + V^PaPp ~ P^PaV^P ~ P^PpVau ~ PvPaVv/3 
- PvPpriaiJ, + P ■ P 1 (riiiaripv + ri^priva - y^riap) • (5.20) 

The first tensor is the result of varying the metric tensor inside the logarithm, while the 
second comes from the metric tensors in the rest of the action. Notice that we enforce both 
transversality and on-shellness except in non-analytic expressions that are infrared singu¬ 
lar. Apart from being gauge-invariant, the above tensors respects local energy-momentum 
conservation 


q IJ V IWQ g = (fV ll:VOl Q = 0 

q^K^P = = 0 • (5.21) 

Indeed this property is guaranteed for the tensor M^ ua p since it is the variation of a local 
operator, but it is gratifying to see that the same applies for T>^ va p which is the variation 
of a purely non-local object. 

6 The 1Z 3 action 

In this section, we perform the matching procedure outlined in section 4. It is more conve¬ 
nient to work in momentum space, and so we list the momentum-space expansions of the 
different curvature invariants in an appendix. 
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6.1 Terms including 1/V 2 


Here we display the non-linear completion of the anomalous contribution to the effective 
action. At the linear level, we had [22] 

T„iM,h\ = I ( ^ 4 ( 6 . 1 ) 

where 

= 19 2 7r 2 (PaPB - P ■ p'Vap)(Q^Qu ~ Q^u + Q 2 p^) ■ (6.2) 

The non-linear completion commences by proposing the ansatz 

( 3 ) p / I -1 

r po i e [g,A] = J d 4 x^ [P S F^F^^R + P Rlc F' 6 ljF a i J ^ 2 R a 0 

+ pRiem F( P F ^^_ R a\ ( 8 . 3 ) 


where the choice of the form factor is easily motivated by the presence of the massless pole 


1 1 

^g 2 ^ □ 


(6.4) 


Using the expansions provided in the appendix, one can form a linear system to solve 
for the three coefficients. It naively appears that the system is overdetermined since the 
expansion of the curvature invariants contain tensor structures that do not appear in eq. 
(6.2). Nevertheless, one only finds exactly three independent equations which uniquely 
yields 


P S 


_ pRic 

192vr 2 ’ 


1 pRiem 

48tt 2 ’ 


1 

96-7T 2 


(6.5) 


We can use the Weyl tensor to change the curvature basis which is very useful to discuss 
the conformal (non)-invariance of the effective action. In 4 D, the Weyl tensor reads 


1 R 

Cfj.ua/3 = ^(^SfiaRufi QfipRva 9uaRfj.fi + 9uftRfia) H jV {.9fia9uj3 9fif)9ua) 


( 6 . 6 ) 


Hence, eq. (6.3) becomes 


( 3 ) /* 11 
I \ole[9, A] = I d A x^-g(P s F, v F^^R + P c F J F^ (6.7) 


where 


pS = 1 pC _ 1 

5767r 2 ’ 967T 2 


V 2 


( 6 . 8 ) 


In fact, the coefficient of the Ricci scalar piece is indeed related to the beta function of the 
theory as could easily be checked by consulting the effective action in fermionic QED [22]. 
One finds 

bi 1 1 


P - 12’ bboson ~ 487T 2 ’ bf 


ermion 


12 vr 2 


(6.9) 
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6.2 Terms including (logV 2 )/V 2 

In the linear action, we also found a logarithmic non-locality which reads [22] 

m.ft] = -\J ( 2^4 log (=^j (6.10) 

where A4j [ iya p has been given previously in eq. (5.20). Although the appearance of Ai° 
might suggest that the above action could be matched onto the quadratic basis, this is in 
fact impossible. We show next that the action can only be matched onto the cubic basis 
with the following form factor 

T log [g, A] = J d\^~g (V ^ R + L Ric R a0 

+ L me m F ^og(^Jg 2 ) Ra ^ ( 6 .H) 


The 1 /□ is inserted for dimensional consistency at this stage as it comprises the only possible 
non-local object one can employ. The matching procedure is the only way to decide on the 
consistency of the ansatz. Once again, using the curvature expansions in the appendix one 
ends up with three independent equations which uniquely fixes the coefficients 

L s = j , L Ric = -b s , L Riem = 0 . (6.12) 

The 1 /q 2 factor which results from inserting the inverse d’ Alembertian cancels out against 
factors of q 2 in the curvarure invariants. Using eq. (3.9), one can rewrite the above action in 
a more transparent form which will prove useful in discussing the conformal (non)-invariance 
of the action 

Tiog[g, A] = b s J d A Xy /g Tg, 1 os( ^Vm 2) ruv (6. 13) 


6.3 Counterterms for the logarithm 

Here we display the counterterms that we need to cancel out the O(h) piece that appears 
in the expansion of the quadratic action eq. (5.17). As we show next, these are third order 
in the curvature. There are two independent tensors in eq. (5.17) which should be matched 
onto two different ansatz. For the tensor Af^ uo ,p, the ansatz is the following 


(3) 

Fct.i[g,A] 


I d 4 x Vs C s F^og{V 2 /g 2 )F^^R + C Ric F^\og{V 2 /g 2 )F a ^R af) 


+ C Riem Fj log (V 2 /g 2 )F^ u ^R 


y2 


(6.14) 


A straightforward matching as before yields 


C s 



£jRic 


b s , C Riem = 0 . 


(6.15) 
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Moving to the tensor 'D pva p, we first notice that in the limit p 2 = p' 2 the non-analytic 
structure becomes 


lnp' 2 — In p 2 
lim -^ y — 

p >2^ p 2 p* — pU 



(6.16) 


which enables us to propose the following ansatz 


(3) 


Tct.2[g,A\ = / drxy/g 


T s F pu ^F^R + T mc F^j F ap R, 


+ t c fJ^f^cx„ 


a (3 


(6.17) 


We choose to work directly in the conformal basis, since it is more convenient. The matching 
yields 


rj-iD _ 


rj-'Ric _ Q rj ~iO _ _ 


C 


12 4 

The same result holds for fermions, substituting bf for b s . 


(6.18) 


7 Remarks on the trace anomaly 


In this section we explore the conformal transformation properties of the different terms 
in the action 6 . We find an interesting dichotomy regarding the terms that give rise to the 
anomaly in response to conformal transformations. This requires a separate treatment of 
scale (global) and Weyl (local) transformations. Since the seminal work of Deser, Duff 
and Isham [37], there has been a consistent effort to understand the precise form of the 
non-local effective action that gives rise to gravitational anomalies. In [26], anomalies 
were geometrically classified to fall into two types. Type A anomalies arise from scale- 
invariant actions, i.e. invariant under a global Weyl rescaling. These are unique and strictly 
proportional to the Euler density of the dimension. On the other hand, type B anomalies 
arise from scale-dependent actions' but the local anomaly itself when denstized is invariant 
under local Weyl tranformations. For example, for a massless minimally coupled scalar in 
2D the anomaly reads 


T v = - R 

11 247T 


(7.1) 


whose density ^fgR is indeed the Euler density in 2D. So this is a type A anomaly, and one 
can check easily that the non-local Polyakov action [38] giving rise to the anomaly is scale- 
invariant. Reigert, following Polyakov, constructed a non-local action in 4D by integrating 
the anomaly [25]. However, the Riegert action was criticized in [26, 39, 41] based on several 
reasons while others [42, 43] argued for its validity. 

The QED trace anomaly falls into type B since its denstized version is indeed (locally) 
conformally invariant, and according to the above classification the generating non-local 

6 See a parallel discussion in [44]. 

' This also means that the action carries an explicit dependence on the renormalization scale /x. 
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action should be scale-dependent. We show below that the two non-local structures present 
in the action are required to generate the correct trace relation whether one performs a 
global or local conformal transformation. Remarkably, the different terms have completely 
different behavior under both types of transformations. In particular, the trace relation is 
generated from the logarithmic non-locality under a scale transformation while the massless 
pole non-locality is responsible for the latter under local ones. 

7.1 Weyl transformations 

Let us commence by considering local transformations. Under an inhnitismal transforma¬ 
tion, we have 


Sc 9fiu — 2a(x)g tllf (7.2) 

which leads to the following transformation of the Christoffel symbol 

= S x V,a + S x V^a - 5 /iI ,V A a- . (7.3) 

From these one readily determines the tranformation of the different curvature tensors. The 
ones we need are 


5 a Rfiu = 2Vfj,V u cr + gf, u V 2 a, 8 a R = 6 V 2 cr - 2aR . (7.4) 

Another object we will need its transformation is the d’ Alembertian operator acting 
on different tensors, in particular, 2 -forms 

5a(V 2 A^) = -2 aV 2 A^ - 2 (V 2 a)A^ - 2(V^a)V x A Xu + 2(V u a)V x A^ x . (7.5) 

where it is understood that A^ is invariant. Once again, let us apply the transformation 
to the quadratic action 

( 2 ) r 

5 CT r[e 2a g,A\ = AL J d 4 x^/g J dm 2 F^{V 2 + m 2 )- 1 (<S ff V 2 )(V 2 + m 2 )" 1 ^ . (7.6) 

Counting the powers of curvature is very important at this stage. The function cr(x) counts 
as a power of the curvature which means that we can freely commute covariant derivatives. 
For example, 

[V M , (V 2 + m 2 )- 1 ] ~ 0{R) . (7.7) 

Using eq. (7.5) and integrating by parts, we find 

S a T[e 2cr g,A] = ^- j d 4 Xy/g j dm 2 F^iy 2 + m 2 ) - 1 (V 2 + m 2 )~ l 

(- 2 a\7 2 F^ - 2 (VV)F + 2aV x V tl F x „ - 2aV x V„F^ . (7.8) 

Now we employ the Bianchi identity 

V^-Fai/ + ^vF^x + = 0 (7.9) 
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to find 


(2) h f 1 

6 a T[e 2 °g,A] = —± J d A x^(y 2 a)F^^F^ 


(7.10) 


Although a prescription to integrate over (dm 2 ) might not seem obvious with the inverse 
operators present in eq. (7.8), one could easily check the above equation by linearlizing 
eq. (7.8) around flat space. It is very important to notice that the above computation 
clearly shows that under the local transformation the log piece does not give rise to the 
anomaly as it does not possess the correct pole structure. Moreover, we show next that eq. 
(7.10) cancels identocally against the contribution coming from the transformation of the 
counter-term. 

Indeed we need not worry about terms containing the Weyl tensor. Moreover, from the 
transformation listed in eq. (7.4) one easily finds 



(7.11) 


given that the field strength is on-shell 


= 0 . 


(7.12) 


The other counter-term transforms as 



(7.13) 


exactly cancelling eq. (7.10) as promised. 


Lastly, the massless pole non-locality of eq. (6.7) is the piece that yields the correct 
trace. To this order in the curvature we only need to keep the 8 a R = 6 V 2 ct + ... term in 
the transformation of eq. (7.4), and neglect the variation of 1/V 2 . Doing this yields 



(7.14) 


which yields the desired trace. In order to see this more simply, and make contact with the 
literature, we can show that all corrections to this result are higher order in the curvature 
by employing the Riegert action [25]. By defining the Paneitz operator [45] 


A 4 = V 2 V 2 + 2 V^R^ - \tf a/ R)V v 

O 


(7.15) 


and 



(7.16) 


where G is the Gauss-Bonnet term 


G = R a/3lS Rap^ - 4 R a0 R a p + R 2 


(7.17) 
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we can see that the Riegert form of this action 

T R [g,A] = J d 4 x^(P s F^F^^-n) (7.18) 

is equivalent to the first term of eq. (6.7) up to terms which are higher order in the curvature, 

(3) 

r#[<7, A] = T po i e [g, A] + 0(F 2 R 2 ). With this form, one can show without approximation 
[45] that 

<5cr-r— = 0 , 5 a 1Z = 6 A 4 (t 

A4 

yielding 

5 a T R [g,A] = ~\j d 4 x^aF^F^ . (7.20) 

The expansion in the curvature has yielded a term which, to this order in the curvature, is 
equivalent to the Reigert action. 

Now we know that a conformal variation of a generic action reads 

<5 CT S = ~ f d 4 X\/g a TJ 1 

and thus indeed eq. (7.20) (likewise eq. (7.14)) yields the correct trace relation. 

7.2 Scale transformations 

A global scale transformation can take a couple of forms. One involves the scaling relations 
shown in eq. (2.10). It is simple to see that this transformation leaves all terms invariant, 
except the covariant logarithm. The logarithmic terms inside the square brackets [...] of 
eq. (2.14) are both shifted by InV 2 —> InV 2 — In A 2 , but In A 2 cancels out leaving the 
whole expression invariant. So in contrast to the above Weyl transformation, this form of 
rescaling yields an anomaly that comes from the covariant logarithm. 

Interestingly in the presence of the metric, there is another way to achieve a global 
scale transformation. In this case the transformation on the metric acts as follows 

9[iv ~> (7.22) 

where a is a constant, not necessarily infinitesimal. This may seem like a subcase of the 
Weyl transformation, but in fact it is distinct [46]. Computationally, a distinction arises 
in that derivatives of a vanish, so that many of the integration-by-parts steps from the 
previous section are not available. 

In this case, the transformation properties of the different curvature tensors proceeds 
easily 


(7.21) 


(7.19) 


R, 


■fiv 


R 


plSl 


R^e-^R, C%-K7'I 


va(3 


(7.23) 
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With these relations in hand, we can apply a scale transformation to the covariant action 
to recover the trace relation. We start with the quadratic action 

( 2 ) h. r 

5 a r[e 2rT 9l A] =-a ^ J d*Xy/g F^F^ . (7.24) 

All terms with the form factor 1/V 2 are scale invariant, hence 

(3) (3) 

5 a T pole [e 2(T g,A} = 0, d^T ct . 2 [e 2a g, A] = 0 (7.25) 

while terms with the form factor In V 2 /V 2 cancel each other identically as described above 

(3) (3) 

S a r log [e 2a g,A] = -S a T ct . 1 [e 2cT g,A} . (7.26) 

The anomalous trace of the energy-momentum tensor is easily determined from eq. (7.21) 
and hence eq. (7.24) correctly reproduces the trace relation 

V = ■ ( 7 - 27 ) 

Again it is the logarithm which is the determining factor for the anomaly. 

8 Summary 

We have used a method which we refer to as non-linear completion in order to match 
the one-loop perturbative expansion of the QED effective action to a covariant expansion 
in the generalized curvatures. Within this procedure, the matching has been unique and 
relatively simple to implement. The results are given in eqs. (2.17), (2.14) and (2.15). 
These summarize the one-loop perturbative calculation involving one gravitational vertex. 

The effective action also encodes the anomaly structure of the theory. For the anomaly, 
the important aspect is to generalize the feature that appeared as lnD in flat space. Our 
generalized result Eq. 2.14 contains many terms when expressed in terms of covariant 
derivatives and curvatures. All of these are required in order to both match the one loop 
perturbative calculation and to respect general covariance. There is also an interplay be¬ 
tween these terms and various forms of scale and/or conformal invariance. There is a 
dispute in the literature about whether the anomaly comes from logarithmic terms or from 
the Riegert action[25], e.g. see[26, 39, 41] and [42, 43]. In our explicit computation, we 
showed that both forms are required in order for the action to respond properly to different 
types of transformation. 

Given the simplicity of the perturbative result eq. (2.7), and the complexity of the 
expansion in the curvature eq. (2.14), one suspects that there is a better covariant repre¬ 
sentation for this result. However, the expansion in the curvature is one of the few covariant 
approximation schemes available and therefore needs to be well explored. We are not pre¬ 
pared to propose an improved representation in this paper, and are only trying to match 
the perturbative result to the standard form found when performing an expansion in the 
curvature. We (hopefully) reserve this improved representation to a future publication. 
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In addition, we note that some of the higher order terms in the curvature expansion 
have the potential to be singular in the infrared, and these higher order terms have only been 
lightly explored. More work is in progress to understand whether the non-local expansion 
in the curvature is useful in phenomenological applications [47]. 
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A Momentum space representation 

In this appendix, we collect all the momentum space representations of the different curva¬ 
ture invariants. The quadratic density ~gF 2 is very simple and we do not list here. Moving 
to the cubic invariants, we find 


/ 


d 4 x^F^F^R 


/ 



where 


= 2 (p • p'va/s - P 0 p' a ){q^qu - q 2r n 


(A.2) 


The invariant including the Ricci tensor reads 



where 


„ = -p ■ v [(Q„Q„ + q^Vafi ~ 2 {p'^PpVav + rfvPpVan + PaP^pv + PaPvVpv) 

- q 2 {jl^rivp + rival !+ VPaVfiv) - IppPaVnv) ~ ^^PpPa • (A.4) 


Lastly, the invariant including the Riemann tensor reads 



where 



(A.6) 





In the above, Q = p + p' and q = p — p'. Moreover, transversality and on-shellness are 
assumed 


p 2 = p ' 2 = 0, p ■ A(p) = p’ • A(p') = 0 . (A.7) 

One can easily check that the matrix elements are both gauge-invariant and satisfiy 
energy-momentum conservation. For example, 

P a M s a ^ u = = 0 . (A.8) 

Moreover, one can use the Weyl tensor given in eq. (6.6) to get 

I ^Xy/gFjF^C^ = J £ ^ )A A a {p)A^-p')h^{-q) + 0(h 2 ) (A.9) 

where 

= g ip- priced - PpPa) ( ( hAv ~ 3 Q^Qu ~ q 2 V^) • (A. 10) 

The above tensor is clearly traceless as required since it stems from a conformally invariant 
Lagrangian. 
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